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Abstract 

Quantum annealing is a proposed combinatorial optimization technique meant to exploit 
quantum mechanical effects such as tunneling and entanglement. Real-world quantum annealing- 
based solvers require a combination of annealing and classical pre- and post-processing; at this 
early stage, little is known about how to partition and optimize the processing. This article 
presents an experimental case study of quantum annealing and some of the factors involved in 
real-world solvers, using a 50A-qubit D-Wave Two machine and the graph isomorphism problem. 

To illustrate the role of classical pre-processing, a compact Hamiltonian is presented that enables 
a reduced Ising model for each problem instance. On random N-vertex graphs, the median num¬ 
ber of variables is reduced from to fewer than N log 2 N and solvable graph sizes increase from 
N = 5 to N = 13. Additionally, a type of classical post-processing error correction is evaluated. 

While the solution times are not competitive with classical approaches to graph isomorphism, the 
enhanced solver ultimately classified correctly every problem that was mapped to the processor 
and demonstrated clear advantages over the baseline approach. The results shed some light on 
the nature of real-world quantum annealing and the associated hybrid classical-quantum solvers. 

Quantum annealing (QA) is a proposed combinatorial optimization technique meant to exploit 
quantnm mechanical effects such as tunneling and entanglemeniP. Machines purportedly imple¬ 
menting a type of qnantum annealing have recently become available. While the extent of qnan- 
tumness in these implementations is not fully understood, some evidence for quantum mechanical 
effects playing a useful role in the processing has been appearin^^HSl Aside from the debate over 
quantnmness, there are interesting questions regarding how to effectively solve a real-world problem 
using a quantum annealer. Quantum annealing-based solvers require a combination of annealing 
and classical pre- and post-processing; at this early stage, little is known abont how to partition 
and optimize the processing. For instance, cnrrent quantum annealers have severe practical limita¬ 
tions on the size of problems that can be handled. Can the pre-processing algorithms be modified 
in order to improve scalability? A second question involves post-processing. Quantum annealers 
provide solutions to an embedded version of a problem involving physical qubits. Post-processing is 
generally required for translating these to solutions to the original problem involving logical qubits 
(aka variables). Occasionally, a chain of physical qubits representing a single variable resolves to 
an inconsistent state, a scenario known as a broken chain. Studies are needed regarding broken 
chains and the possibility of classical error correction dnring post-processing. 
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This article presents an experimental case study of quantum annealing and some of the factors 
involved in real-world solvers, using a 504-qubit D-Wave Two machine. An example of parsimonious 
pre-processing is considered, along with post-processing error correction. Through experiments on 
a 504-qubit D-Wave Two machine, we quantify the QA success probabilities and the impact of 
the methods under study. We use the graph isomorphism (GI) problem as the problem of focus. 
The GI problem is to determine whether two input graphs Gi ^2 are in essence the same, such that 
the adjacency matrices can be made identical with a relabeling of vertices. This problem is an 
interesting candidate for several reasons. First, an accurate quantum annealing-based solver for 
GI has never been implemented. Second, quantum approaches can sometimes provide new insight 
into the structure of a problem, even if no speedup over classical approaches is achieved or even 
expected. Third, the GI problem is mathematically interesting; though many sub-classes of the 
problem can be solved in polynomial time by specialized classical solvers, the run time of the 
best general solution is exponential and has remained at since The classical 

computational complexity of the problem is currently considered to be NP-intermediat^, and the 
quantum computational complexity of the problem is unknown. Graph isomorphism is a non- 
abelian hidden subgroup problem and is not known to be easy in the quantum regimdSEH Lastly, 
the GI problem is of practical interest. It appears in fields such as very large-scale integrated circuit 
design, where a circuit’s layout graph must be verified to be equivalent to its schematic grapli^, 
and in drug discovery and bio-informatics, where a graph representing a molecular compound must 
be compared to an entire database, often via a GI tool that performs canonical labeling. 

This article relates to previous works as follows. A pre-print by King and McGeoch discusses 
tuning of quantum annealing algorithms, including the use of low-cost classical post-processing 
error correction similar to what is evaluated in this articl^^. Our study goes further regarding 
pre-processing (designing a Hamiltonian to generate compact Ising models) and covers graph iso¬ 
morphism rather than problems such as random not-all-equal 3-SAT. A work by Rieffel et al. maps 
real-world problems such as graph coloring to a D-Wave quantum annealei^. Regarding the graph 
isomorphism problem in particular, multiple attempts have been made using adiabatic quantum 
annealing. One of the first attempts assigned a Hamiltonian to each graph and conjectured that 
measurements taken during each adiabatic evolution could be used to distinguish non-isomorphic 
pair^. A subsequent experimental study using a D-Wave quantum annealer found that using 
quantum spectra in this manner was not sufficient to distinguish non-isomorphic pair^^^. A second 
approach converted a GI problem to a combinatorial optimization problem whose non-negative cost 
function has a minimum of zero only for an isomorphic pair. The approach required AI[log 2 iV] 
problem variables and additional ancillary variables. It was numerically simulated up to A = 7 but 
not validated on a quantum annealing processoi^. An alternative GI Hamiltonian was proposed 
by Luca^. 

Proposed Solver 

Preliminaries. Before proceeding, we briefly define some key concepts. In quantum annealing, a 
system is first initialized to a superposition of all possible states as dictated by Hamiltonian Hinu, 
then slowly evolved such that the initializing function becomes weaker and the energy function 
of interest, defined by the problem Hamiltonian Hp, becomes dominant. The time-dependent 
combination of the two Hamiltonians is referred to as the system Hamiltonian H 


H (t) — A (t) Hinit -|- B (t) Hp, t G [0, T], 


( 1 ) 
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where A (t) and B (t) are monotonic functions representing the annealing schedule up to time T. 
In this work, we focus on problem Hamiltonians Hp. The problem Hamiltonian represents an Ising 
model, where spin variable Si G {—1,1} are subject to local fields hi and pairwise interactions with 
coupling strengths Jij^- 


Hp — ^ ^ hiSi + ^ ^ JijSiSj. (2) 

i i,j 

It is often the case that an Ising problem [h, J] involves variable interactions not supported 
by a quantum annealing processor graph. An example is a variable with degree higher than the 
maximum supported by the D-Wave Chimera architecture (six). One common strategy is to find a 
minor embedding of the problem graph in the processor graphP^, in which a set of physical qubits 
is used to represent each variable. Each qubit is strongly coupled to at least one other qubit in the 
set, in an effort to keep the entire set in a consistent state. These sets are commonly referred to as 
chains. We refer to a minor-embedded version of a problem as an embedded Ising problem [h', J']. 

Baseline Hamiltonian. We first describe a baseline penalty Hamiltonian for the GI problem, 
building upon the Hamiltonian described in Luca^^. The problem input is a pair of simple, 
undirected Wvertex graphs Gi, 2 . Penalties are applied such that the ground state energy is zero 
if the pair is isomorphic and greater than zero otherwise. The intent is for an energy minimization 
process (such as quantum annealing) to provide a solution to this decision problem. 

A binary variable Xu,i is created for every possible mapping of a vertex u in G 2 to a vertex 
i in Gi; in a solution, the variable is 1 if u is mapped to i, and 0 otherwise. Since the model 
is restricted to pairwise interactions over binary variables, it represents a quadratic unconstrained 
binary optimization (QUBO) formulation, which can be readily converted to an Ising model. There 
are two types of penalties that can be applied to an interaction between two variables. One type 
(Cl) penalizes vertex set mappings that are not bijective, for instance a mapping in which both xiq 
and X 2 ,i are set to 1 (since vertices 1 and 2 in G 2 cannot both map to vertex 1 in Gi). The term in 
the QUBO matrix corresponding to the xiqX 2 ,i interaction is thus set to Gi. Note that as long as 
at least one of the two variables resolves to 0, the energy of the coupling (Gixi, 1 X 2 , 1 ) is zero. The 
second penalty (G 2 ) applies to edge inconsistencies. An example of an edge inconsistency is when 
vertices u,v are mapped to i,j, and uv G E 2 while ij 0 Ei. Additional details regarding this style 
of Hamiltonian can be found in Luca^^. 

In Luca^i^, couplings can incur either zero, one, or two penalties. Here, we require that couplings 
be penalized no more than once, in order to achieve a simple set of coupler values (e.g. { 0 , 1 } 
instead of {0,1,2}) amenable to quantum annealing. This is enforced by setting Gi = G 2 > 0, 
and preventing double penalties. An edge-related penalty is applied to a coupling only if there is 
not a vertex mapping penalty. As an example, if vertices u,v are mapped to i,j and i = j, then 
the bijection has been violated and coupling Xu,iXv,j will incur a vertex mapping penalty (Gi); 
therefore, an additional edge inconsistency penalty (G 2 ) need not be applied. Edge inconsistency 
penalties are only applied if f 7 ^ j and u ^ v. The complete baseline Hamiltonian is 





+ C 2 ^ 




^ ^ Xu,iXv,j T G 2 ^ ^ 

u,v&E 2 i,j&Ei 



U,V^E2 


( 3 ) 
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where the i ^ j and u ^ v conditions represent the main modifications to Luca^^. This QUBO 
form can be converted to an Ising form using the relation Xu,i = • 

Compact Hamiltonian. The approach embodied in the baseline Hamiltonian Hi suffers from 
a severe lack of scalability. For A^-vertex input graphs, it requires N'^ logical variables. Moreover, 
due to the limited direct connections between qubits in the D-Wave Chimera architecture, problems 
are often given a minor embedding into the processor working graph. This typically involves repli¬ 
cating variables across multiple qubits. Thus, the qubit requirements can reach O (fV^). Problems 
mapped in this way to a ~ 500-qubit processor tend to be limited to N = 5 or 6. We now investi¬ 
gate whether a more effective Hamiltonian can be designed. The idea is that many variables and 
interactions are unnecessary, and information indicating so can be leveraged up front during the 
requisite pre-processing. Note that an isomorphic mapping requires the vertices in each matched 
pair to have the same degree. Thus, degree information can be used to decide whether two vertices 
are eligible to be matched. We propose a compact Hamiltonian H 2 that avoids creating variables for 
vertices of different degree. A second, minor simplification deals with isolated vertices {degree = 0). 
If Gi ^2 each have k isolated vertices, an isomorphic mapping of such vertices is trivial and thus no 
variables or penalties for those vertices need be modeled. If Gi ^2 have a different number of isolated 
vertices, then they also have a different number of non-isolated vertices and existing variables and 
penalties for those will suffice. Thus we only create variables and penalties for vertices with degree 
greater than zero. 

Given the two enhancements, the total number of variables required is where \di\ is 

the multiplicity of the set di containing vertices of degree i. In the worst case of regular graphs of 
degree r greater than zero, all nodes have the same non-zero degree and thus the simplifications 
provide no benefit— N'^ variables are still required. Many real-world graphs are not regular and for 
these the benefits can be large. An example illustrating the concept is provided in Fig. 
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Figure 1 |Example of variable reduction. 
The entire compact Hamiltonian H 2 is 
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V {deg{i) = deg{u) > 0 ; deg{v) = deg{j) > 0 ) : 

u \ i / i \ u / 


+ E 


^ ^ Xu,iXv,j 
u,v&E2 


+ E 


^ ^ ^u,i^v,j ’ 
U,V^E2 

U^V 


H2 was validated using a software solver (D-Wave’s ising-tieuristic version 1 . 5 . 2 ) that provides 
exact results for problems with low tree width. In exhaustive testing of all 2 ^^ N = 4 pairs and 
2^0 AT = 5 pairs, the ground state energy of H2 was confirmed to be zero for isomorphic cases and 
greater than zero for non-isomorphic. 

As a part of a quantum annealing-based solver, an algorithm can be employed that accepts 
a graph pair as input and uses the Hamiltonian to generate an associated QUBO problem (later 
converted to an Ising problem). The proposed algorithm using the compact Hamiltonian H2 is 
presented (in pseudo code) in Algorithm 


Input: adjacency matrices Gi^2 
Output: QUBO matrix Q 

// Define QUBO variables 

for each node v in G2 do 
for each node i in Gi do 

if deg{v) == deg{i) and deg{i) > 0 then 
I create QUBO variable 
end 
end 
end 

// Populate QUBO matrix with penalty terms 

for each pair of different QUBO variables Xu,i and Xvj do 

// Assign penalty for node mapping conflict 

if i == j then 

I penalize mapping U —>■ i, V — ^ j II Two nodes in G2 map to the same node in Gi 

end 

else if u == v then 

I penalize mapping U ^ i, V ^ j II A node in G2 maps to two nodes in Gi 

end 

// Assign penalty for edge discrepancy 

else if Gi{i,j) ■= G 2 {u,v) then 

I penalize mapping U —>■ i, V —>■ j II Edge in one graph, non-edge in other 

end 

end 

for each QUBO variable do 
I Assign value along diagonal of Q 

end 

Algorithm 1: High-level algorithm for generating a QUBO problem using the 
proposed compact Hamiltonian H2. 

Complete Solver Flow. The high-level flow of the proposed GI solver is shown in Fig. The 
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problem input is a graph pair Gi^2- In this article, the graph types considered are random, simple, 
undirected graphs. Graphs are generated using the Erdos-Renyi modeP^ G{N,p) where we set 
the probability of an edge being present p = 0 . 5 . An advantage of G{N, 0 . 5 ) graphs for an initial 
study is that all graphs are equally likely. This type has been used in classical graph isomorphism 
work as welP. In step 1, the input graphs and the Hamiltonian formulation of interest (e.g. Hi or 
H2) are used to generate a QUBO problem which is then converted to an Ising problem [h,J]. An 
example of an algorithm for generating the QUBO problem is shown in Algorithm The Ising 
problem is then compiled to a specific quantum annealing processor in step 2 . A main task is to 
find sets of physical qubits to represent the problem variables (aka logical qubits); this is achieved 
by providing the J matrix and the processor working graph to the D-Wave f indEmbedding () 
heuristic!^. Subsequently, the parameters of the embedded Ising problem [h', J'] are set following 
certain strategies such as the use of a random spin gauge (see Methods). The embedded Ising 
problem, sometimes referred to as a machine instruction, is submitted to the quantum annealing 
machine along with several job parameters. The quantum annealing job is executed in step 3 
and solutions are returned in the form of strings of two-valued variables. These solutions and 
energies are associated with the embedded problem, not the original Ising problem. Therefore a 
post-processing step is necessary (step 4 ), in which the state of each qubit chain is plugged into 
the cost function of the Ising problem. A difficulty arises when the states of the qubits in a chain 
are inconsistent, a case referred to as a broken chain. In the proposed solver, broken chains can be 
handled by either discarding the associated solution, or by performing majority voting over each 
chain. The two strategies are compared empirically in Results. Given a solution to the original 
Ising problem, the solution energy can be calculated. If the lowest energy is zero, then the input 
pair can be declared isomorphic and no further jobs are necessary. Otherwise, a decision must be 
made whether to repeat the process from step 2 or to stop and declare that isomorphism could not 
be established. 
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Problem input: Designed Hamiltonian 

graphs 2 (e.g. H 2 ) 



Figure 2 | Graph isomorphism solver flow. 


Results 

Ising Model Scaling. To compare the resource requirements of the two proposed Hamiltonians, 
100 pairs of graphs are used as inputs to Step 1 of the solver flow (Fig. [^, where 50 pairs are 
isomorphic and 50 are non-isomorphic for each size up to = 100. Since Hi models a variable for 
each possible vertex pair, iV^ variables are required by definition. Ising problems generated using 
H 2 are found to use fewer variables than Hi] scaling of the median problems fits to 0.748Ai^^'^®. 
Incidentally, this indicates that most problems have fewer variables than with the Gaitan et al. 
approach, which entails iV[log 2 iV] plus ancillary variabled^. The variable scaling is illustrated in 
Fig. i In addition to the number of variables, a second resource metric is the number of non¬ 
zero interactions between variables; dense interactions make the minor embedding problem more 
difficult. We find that the scaling of variable interactions has been improved from O for Hi 
to O for H 2 (where B? = 0.9991). 

Embeddability. Next, we compare the embeddability of the two approaches, in other words 
the extent to which Ising problems can be minor-embedded in a given processor graph. The 
processor of choice is the D-Wave Two Vesuvius-6 processor housed at USC ISI. At the time of 
this writing, the working graph contains 504 qubits and 1427 couplers. Embedding is attempted 
using the D-Wave f indEmbeddingO heuristicP^ with default parameter values such as 10 tries per 
function call. As shown in Fig. embeddings are found for the majority of problems only for 
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sizes N < 6 when using Hi, but sizes N < 14 with H 2 (Fig. |^). The median number of qubits 
across all problems scales as O for Hi and has been reduced to O (iV^'^®) for H 2 (Fig. |^)- 



Figure 3 | Scaling of the number of Ising model variables. Hi requires variables by construction; 
H 2 scales more efficiently—the median variable requirement fits to y = 0.748N^ '^^{R^ = 0.9995). Bars 
indicate the maximum and minimum. For reference, the green dotted line represents A^|"log 2 N~\. Inputs are 
50 isomorphic and 50 non-isomorphic pairs of G{N, 0.5) graphs. 
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Figure 4 jEmbeddability when targeting the USC-LM Vesuvius processor’s 504-qubit, 1427- 
coupler working graph. The inputs are 50 isomorphic and 50 non-isomorphic G{N, 0.5) pairs at each 
size, (a) Fraction of problems for which an embedding was found within 10 tries, using the D-Wave 
f indEmbeddingO heuristic, (b) Median number of utilized qubits among all 100 problems. Hi data points 
fit to 2 / = 0.235iV4'22 (dotted line; = 0.9996); H 2 fits to 0.0723iV3-29(i?2 ^ 0.9955). The 65'^'' and 35“^ 
percentiles are marked with bars (not always visible). 













































Experimental Quantum Annealing for Graph Isomorphism. The accuracy of the solver 
described in the previous section was measured via trials conducted on a D-Wave Two Vesuvius 
quantum annealing processor. Several alternative strategies were compared—the use of Hamiltoni¬ 
ans Hi vs. H 2 , running a single job per problem vs. multiple jobs, and the use of error correction. 
Note that by construction of the Ising models using a penalty Hamiltonian, problems with non¬ 
isomorphic input graphs cannot achieve a zero energy state, regardless of annealing results. The 
main challenge for the solver is to find the zero energy state for isomorphic pairs. Thus, we first 
focus on the isomorphic case. One hundred isomorphic pairs were input into the solver for each 
size N from 3 to 20. 

For one strategy in particular the zero energy state was always eventually achieved—the use of 
Hamiltonian H 2 combined with multiple jobs and error correction. Thus, with this strategy there 
were no false negatives and classification accuracy reached 100% of the embeddable problems, as 
shown in Table For the most difficult problem, the zero energy state was achieved on the 9**^ 
job. All other strategies incurred false negatives. For the successful strategy, the expected total 
annealing time was calculated (as described in Methods). Results are shown in Fig. 

Table 1 | Number of isomorphic-input problems embedded and correctly classified as isomorphic 
via quantum annealing. One hundred problems were attempted at each problem size. All embedded 
problems were solved when using H 2 , error correction, and multiple jobs. 



Size of input graphs (number of vertices) 


3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

Number of problems for 
which an embedding was 
found using H\ 

100 

100 

100 

74 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Number of problems for 
which an embedding was 
found using H 2 

100 

100 

100 

99 

99 

96 

95 

93 

87 

83 

54 

39 

20 

7 

4 

1 

Number of problems 
solved using Hi, no 
error correction 

100 

99 

98 

13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Number of problems 
solved using Hi, error 
correction 

100 

100 

99 

43 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Number of problems 
solved using H 2 , no 
error correction 

100 

100 

100 

98 

96 

93 

86 

75 

67 

53 

34 

23 

9 

1 

2 

0 

Number of problems 
solved using H 2 , error 
correction 

100 

100 

100 

99 

99 

96 

90 

80 

73 

65 

40 

25 

12 

3 

2 

1 

Number of problems 
solved using H 2 , error 
correction, multiple jobs 

100 

100 

100 

99 

99 

96 

95 

93 

87 

83 

54 

39 

20 

7 

4 
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Figure 5 |Total expected annealing time when using Hamiltonian H 2 , multiple jobs, and classi¬ 
cal majority voting. Based on the quantum annealing results for 100 isomorphic-input problems. Missing 
data points indicate that no embedding was found and quantum annealing was not attempted. 


For completeness, non-isomorphic pairs were run as well, using Hamiltonian H 2 and error cor¬ 
rection. Since in the worst case nine jobs were required to correctly classify the isomorphic pairs 
above, nine jobs were submitted for each non-isomorphic problem. One hundred non-isomorphic 
G{N, 0.5) problems were tested at each size between N = 3 to 14; of the 1200 problems, embed¬ 
dings were found for 1186. In addition, pairs of isospectral non-isomorphic graphs (PINGs) were 
tested. All N = 5 PINGs were tested (150 permutations), as well as 100 random N = 6 PINGs. 
As expected, none of the non-isomorphic problems achieved a zero energy state and thus none were 
classified as isomorphic. In other words, there were no false positives. 


Discussion 

Several observations can be made from this case study. First, the formulation of the cost func¬ 
tion (Hamiltonian) can have a noticeable impact on quantum annealing results. For the graph 
isomorphism problem, the baseline approach (embodied in Hamiltonian Hi and in Luca^^ blindly 
creates QUBO variables for every possible vertex pair, whereas the proposed Hamiltonian H 2 is 
more parsimonious. Variable requirements reduced from to fewer than AIlog 2 A^ (Fig- on 
the graph type under study, allowing larger problems to be solved (Fig. and Table [^. Along 
with RieffeP^, this is one of the first quantum annealing studies to experimentally quantify the 
effect of alternative Hamiltonian formulations. One of the impacts of this observation is increased 
appreciation for the fact that all quantum annealing-based solvers are actually classical-quantum 
hybrids and that focus must be placed on effectively partitioning the processing and optimizing the 
classical portion. A caveat is in order—if the classical side is made to do too much work then the 
quantum annealing aspect becomes trivial and of little value. Further work is needed in identifying 
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the specific strengths of annealing processors, and in leveraging the two sides appropriately. 

A second observation is that using majority voting error correction during post-processing can 
in some cases provide a benefit. Previously, such majority voting was evaluated for a different 
set of problems (scheduling) and was not found to provide a significant benefilP^. In our context, 
there were many problems for which the zero energy ground state solution was only achieved 
when using this post-processing; without this form of error correction (in other words, when all 
solutions containing a broken chain were discarded), false negatives occurred. For instance, at 
N = 12, 53 of 83 embedded problems were solved on the first job without using error correction, 
and an additional 12 problems were solved by applying error correction (Table [^. Classical error 
correction strategies other than majority voting should be explored and assessed in future studies, 
and their costs quantified. 

To our knowledge, the evaluated solver is the first validated, experimental implementation of 
a QA-based graph isomorphism solver. While it ultimately classihed every embeddable problem 
correctly and demonstrated clear advantages over the baseline approach, it has serious limitations as 
a graph isomorphism solver. The problem sizes are not competitive with those handled by classical 
solvers, which can handle G {N, 0.5) graphs with thousands of vertice^ and even for the hardest 
graph types can handle hundreds of vertices before running into difficult}^^. Similarly, the scaling 
of the total annealing times (Fig. is not competitive with classical scalin^^. Ultimately, new 
approaches are likely needed if quantum annealing is to contribute to graph isomorphism theory 
or practice. Fortunately, the case study provides some new insight into experimental quantum 
annealing, and contributes methods that have relevance beyond the GI problem. It is hoped that 
the experimental evaluation of alternative Hamiltonian formulations adds to the understanding 
of the factors affecting quantum annealing performance, and that the demonstration of majority 
voting raises new questions about the role of post-processing for a variety of problems. 


Methods 

Quantum annealing experiments were performed on the D-Wave Two machine housed at USC ISI and op¬ 
erated by the USC-Lockheed Martin Quantum Computing Center. Experiments were conducted in October 
and November, 2014. The working graph of the machine’s Vesuvius-6 quantum annealing processor consisted 
of 504 qubits and 1427 couplers during this period. The pattern of working qubits is shown in Figure The 
qubit temperature was estimated to be 16 ± 1 mK. Additional processor specifications include a maximum 
anti-ferromagnetic mutual inductance of 1.33, and j amplitude of 7.5 ± 1 
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Figure 6 | Physical layout of the working qubits in the USC-LM D-Wave Two Vesuvius-6 
processor as of October 10, 2014. 504 working qubits (gray), 8 non-working (black). All 1427 couplers 
connecting working qubits are part of the working graph. 

Simple undirected V-vertex graphs were constructed according to the Erdos-Renyi G (ri,p) modeP^with 
n = N and with the probability p of including each edge equal to 0.5. Non-isomorphic pairs were generated by 
creating two graphs as above and checking for non-isomorphism using the MATLAB graph! somorphismO 
function. Isomorphic pairs were generated by generating a single graph then applying a random permutation 
to arrive at the second graph. For each pair of input graphs, an Ising model was created using equation 
( [^ or (|^. Programming was performed using MATLAB R2014a win64 and the D-Wave MATLAB pack 
1.5.2-beta2. The current version of the D-Wave sapiFindEmbeddingO function cannot embed Ising models 
with more than one connected component (i.e. a set of variables that interact only with each other and not 
any of the remaining variables); therefore, models with this characteristic were not included in the input 
data. When attempting to generate 100 input pairs for each size, such disconnected models occurred no 
more than 4 times for each size N > 14. Similarly, the heuristic cannot accept models with fewer than two 
variables, so in the rare case of a trivial Ising problem with fewer than two variables (e.g. a non-isomorphic 
pair with no matching degrees), dummy variables were added to the problem. 

The hi values of the Ising problem were split evenly across each qubit in the associated chain in the 
embedded Ising problem. The Jy values of the Ising problem were assigned to a single coupler connecting 
two variable chains in the embedded problem. The magnitudes of the embedded h' and J' were scaled 
together such that the maximum magnitude reached 20% of the full range supported by the processor; the 
range of the embedded h[ values was [—0.4, 0.4] and the range of the embedded Jb values coupling different 
variables was [—0.2, 0.2]. This 20% value was determined empirically to provide good performance on the 
median difficulty problem at the largest sizes. Subsequently, the Jb values connecting physical qubits within 
a chain were set to the maximum ferromagnetic value (—1). A single random spin gauge transformatiorP was 
then applied to each embedded problem, with a gauge factor G {—1,1} associated with each qubit and 
transformation h'i -G aih'^^Jh —)■ aiUjJh. One job was submitted to the quantum annealer per embedded 
problem; some Ising problems were associated with multiple embedded problems and jobs. After each 
programming cycle, the processor was allowed to thermalize for 10 ms (the maximum supported by the 
machine). The annealing time was set to the minimum value of 20 fis. The number of annealing and readout 
cycles per programming cycle was 40000, which allowed the total job time to be within the limits of the 
machine (Is). The readout thermalisation time was set to the default value of 0. Regarding error correction 
through majority voting of chains of physical qubits, ties were broken by choosing the spin up state. The 
probability of achieving the zero energy state on job k is denoted 
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(5) 


Po,k = 


number of annealing cycles achieving zero energy 
number of annealing cycles 


When multiple jobs are required, we calculate the geometric mean in the style of Boixo et 


al.l^ 


K 

Po = l-l[{l-Po,k)^ ■ ( 6 ) 

k=l 

The total annealing time required to reach 0.99 probability of success was calculated by multiplying the 
annealing time by the expected number of annealing cycles (repetitions R) using the formuloP. 


ln(l - 0.99) 
In (1 - Po) 


(7) 
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